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Abstract. 



The geometrodynamics of the spherical gravity with a selfgravitating thin 
dust shell as a source is constructed. The shell Hamiltonian constraint is de- 
rived and the corresponding Schroedinger equation is obtained. This equa- 
tion appeared to be a finite differences equation. Its solutions are required 
to be analytic functions on the relevant Riemannian surface. The method of 
finding discrete spectra is suggested based on the analytic properties of the 
solutions. The large black hole approximation is considered and the discrete 
spectra for bound states of quantum black holes and wormholes are found. 
They depend on two quantum numbers and are, in fact, quasicontinuous. 
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1 Introduction 



In the absence of the theory of quantum gravity we have to construct a 
new theory every time we want to quantize some classical gravitating object. 
The most interesting models of the kind are black holes and cosmological 
models . Here we are interested in the quantum black holes models. These 
models deserve consideration for many reasons. The black hole physics gives 
us an example of the strong gravitational fields. The existence of the event 
(apparent) horizons causes the Hawking's evaporation of the black holes . 
The fate of the evaporating black holes becomes a subject of interest. The 
quantum theory may throw some light to many problems of the classical 
black hole physics. 

There were many many attempts to construct such quantum models. The 
most interesting of them are described in 0, @].In these papers the quan- 
tum theory of the eternal Schwarzschild black hole was constructed . The 
authors introduced many useful and important mathematical tools (in the 
present paper the canonical transformation found by Kuchar is widely used, 
but the physical result is rather trivial and obvious. Namely the quantum 
functional depends only on the Schwarzschild masses. The reason for this is 
that the eternal Schwarzschild black hole has no dynamical degrees of free- 
dom. That is , all the matter collapsed classically and all possibly dynamical 
degrees of freedom died in the singularity. In the present paper we treat (or 
try to ) the simplest quantum black hole model. "The simplest" means that 
we consider a spherically symmetric gravity with a self-gravitating thin dust 
shell as a source. We constructed the classical geometrodynamics for the 
system. Quantization of such a model leads to the Schroedinger equation in 
finite differences in the coordinate representation. The shift in the argument 
is along imaginary axes which has very important consequences. One of them 
is that the wave function which are the solution to such an equation should 
be analytical function on the appropriate Riemannian surface. 

It should be noted that it is not the first time we are dialing with the 
finite differences equation . In the toy model constructed by one of us the 
Schroedinger equation in finite differences emerged as a result of the use of the 
proper time quantization. Unlike this toy model the present consideration 
deals with the canonical formalism from the very beginning. Thus , the 
results do not depend on the choice of time. And the appearance of the 
finite differences equation is due to the nonlocal nature of the corresponding 
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Hamiltonian operator. 

In the ordinary quantum mechanics we are deahng with the second order 
differential equations. And we demand that the solution should be at least 
two times differentiable. To find eigenfunctions and spectrum wc need to 
specify a class of functions, usually by imposing appropriate boundary con- 
ditions. In our case of the finite differences operator wc must specify a class 
of function by demanding analyticity (except in the branching points). Our 
experience with the toy model shows that the boundary conditions help us 
to select the wave eigenfunctions (though up to the infinite degeneracy) but 
the are useless in finding the mass spectrum. To find the spectrum we need 
to know only the analytic properties of the solutions, namely their branching 
points. 

It makes our hfe easier. We do not need to solve the Schroedinger equa- 
tion. We should only to investigate the behavior of the solutions in the 
singular points of the corresponding equation. 

The plan of the paper is the following. In the Section 2 wc remind some 
facts from the classical dynamics of the thin dust shells. The classical ge- 
ometrodynamics is developed in Section 3. Section 4 is devoted to the deriva- 
tion of the shell's Hamiltonian constraints. The quantum geometrodynamics 
of the spherical gravity with the thin shell is considered in Section 5. Section 
6 is dealing with the quasiclassical limit of our Schroedinger equation which 
in our case is the same as the large black holes regime. In this section we 
found the quantum black hole and wormhole discrete mass spectra. 
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2 Preliminaries. 



The aim of this paper is to consider a geometrodynamics , both classical and 
quantum of the spherically symmetric gravitational field with self-gravitating 
dust thin shell as a source. 

We start with the description of the model. This is just a self-gravitating 
spherically symmetric dust thin shell, endowed with a bare mass M. The 
whole space-time is divided into three different regions: the inner part ( I ), 
the outer part ( II ) containing no matter fields separated by thin layer III, 
containing the dust matter of the shell. 

The general metric of a spherically symmetric spacetime has the form: 

ds^ = -N^dt^ + L\dr + N'dtf + R^{de^ + sin^ dd(j?) (1) 

where {t,r,6,(l)) are space-time coordinates, N, , L, R are some functions 
of t and r only. Trajectory of the thin shell is some 3-dimensional surface 
S in space-time given by some function f{t): = {{t,r,6,(f)) : r = r{t)}. 
In region I r < f — e, in a region 11 r > f + e, region III is a thin layer 
r — e<r<f + e. 

We require that metric coefficients A^, N^,L and R are continuous func- 
tions but jump discontinuities could appear in their derivatives at the points 
of S when the limit e ^ is taken. 

Contrast to the flat space-time the normal vector to the surface R = const 
may not be only spacelike but also timelike. In the first case the invariant 

F = g''^R,aR,i3 > (2) 

and the corresponding region is called i?-region (fig. 1) (here Qajj is a metric 
tensor, its inverse, R^a denotes the partial derivative with respect to 

the corresponding coordinate, Greek indices run from to 3). In the fiat case 
i?-region occupies the whole space-time. In the second case 

F < (3) 

Such a region is called the T-region (the notions of R- and T- regions were 
introduced in |]T|). It is easy to show that the condition R = (dot denotes 
time derivative) cannot be satisfied in a T-region, hence it should be either 
R > (this region of inevitable expansion is called T+-region), or i? < 
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(inevitable contraction, a T_-region). Correspondingly, it is impossible to 
get R' = (prime denotes the spatial derivative) in i?-regions, and a region 
with i?' > is called an i?_|_-region, while that with i?' < is an R_ region. 
The and -R_- regions correspond to different sides of the Einstein- Rosen 
bridge (see fig. 1). 

The solution of Einstein equations representing the Schwarzschild (spher- 
ically symmetric) black hole is well known and can be put in the form 

ds"^ = -FdT^ + F-^dR^ + R^idO"^ + sin^ ^rf^^) (4) 

where 

f^l-Hf (5) 

and m is the total mass (energy) of the system, G is the gravitational constant 
(equal to the inverse square of Planckian mass, G = MpJ^,in the chosen units 
with c = h = 1; note also that in these units and the radius has dimension 
of inverse mass). 

The metric of 3-dimensional surface , representing the evolution of the 
thin shell can be written as 



ds' 



-dr^ + R{t) {de^ + sin^ edcj)^) (6) 



Here r is the proper time of the observer sitting on the shell. From Einstein 
equations one obtains the equations of motion of the shell in the form 



, - 2 / X 2 QM 

Cin]/ R +Fin — CroutyR + Font = „ (7) 

K 

( see [0]). The quantity a has the following meaning, cr = +1 if radii R 
increase in the outward normal direction to the shell, and a = — 1 if radii 
decrease. Therefore in the i?-regions a does not change its sign, the latter 
being the sign of an i?-region. It is clear now that on "our" side of Einstein- 
Rosen bridge we have a = +1, this we shall call the "black hole case" , 
while on the "other" side a = — 1, this we shall call the "wormhole case". 
On fig. 2-5 we show all possible junctions of inner and outer regions with 
Schwarzschild masses m^^ and rriout with thin shell of bare mass M. This 
figures are rather schematic. The space-time inside the shell (to the left of 
the trajectories ) depends on whether there are other shells inside the given 



6 



one or not. To avoid such a concretization we draw the part of a complete 
inner Schwarzschild space-time with the left infinity. What important here 
are the junction of the space-times and the and the positions of the event 
horizons. 
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3 Canonical formalism for spherically sym- 
metric gravity with thin shell. 

The action functional for the system of spherically symmetric gravitational 
field and the thin shell is 

S = Sgr + Sshell = TfTr / Ry^^d^x + (surfacc terms) — 

" wi+iii (8) 

s 

It consists of the standard Einstein-Hilbert action for the gravitational field 
and matter part of the action describes a thin shell of dust. The surface 
terms in the gravitational action and the falloff behavior of the metric and its 
derivatives were studied in details in |0] . So we will not consider this question 
and will use the results of Kuchar when needed. We will be interested in the 
behavior of the action and constraints on the surface representing the 
shell's trajectory. 

The complete set of degrees of freedom of our system consists of the set 
of N{r,t), N^{r,t), L{r,t), R{r,t) which describe gravitational field and f(t) 
which describes the motion of the shell. 

The metric ([^) has the standard ADM form for 3+1 decomposition of a 
space-time with lapse function , shift vector A^* = (A^^',0,0) and space 
metric hik = diagi^L"^ , E? , siv? 6) given foliation of the manifold on space 
and time. The scalar curvature density has the form 

N\fh(^^^^R+{{TrKf -TrK^))- (g) 
-2 (v^/^),o) + 2 (VhKN' - Vhh'mj) . 

where R and K^^ are the scalar curvature of a space metric hij and exterior 
curvature tensor of a surface t = const. Substituting expression (0) for the 
metric into {^) we obtain the expression for internal and external curvatures 
of the surface t = const in the form 
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and 

Ki=dtag{K:,KlKl) 

k; = ^{l- l'n^ - L{N^y) , (n) 

= < = TO - ^'^') 
Here dot and prime denote differentiation in t and r respectively. 

Contributions to the gravitational action from the terms containing total 
derivatives in (|^) give rise to the surface terms which cancel each other at 
the common boundaries of regions I, II and II, III. So we are left with the 
surface terms at infinity which were extensively discussed in 0. We will turn 
to them later. 

The essential part of the action for gravitational field is just the ADM 
part of the action @ with Lagrangian 

Lgr = Y^iVLi?2 ( (3)i? _ (rrK)2 - TrK') (12) 

Contribution to the action from the integral over the region III in the 
limit e — > is only due to the term containing second derivative of i?,namely 

1 .^.^2^^^ f NRR" f NRIR'] 

——NLR^ ^^^R = - / — — = - / M 13 

IQttG J GL J GL 

III II s 

We will denote by hats variables on S and by [A] = lime^o {A{r + e) — A{f — e)) 
a jump of variable A{r) on the shell surface. 

Substituting the expression (|l|) into the shell part of the action we have: 



Ssheii = -Mj^m-L^ [n- + f) ^dt (14) 
s 

The explicit form of the action (P) with metric (|ID becomes 

i+ii+iii ^ ^ ^ 
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N(k[R']/{GL) + ^m^ + Ti^/L^j- (17) 
N"- (-L [Pl] - Ti) dt 



(16) 



The canonical formalism for this action can be described in the following 
way. Momenta conjugate to corresponding dynamical variables are 

= 6S/6N = 0; 

P^,. = 6S/6N'- = 

Pl = 5S /5L = -^(r'N' -R) 

Pr = 5S /5R=ji^ \R'N'' -R) + -^({LN^y -L 

p^ = 6S/6R = 

Pl = 6S/6L = 

ir = 6S/6r= ^^^^N^ + r) 

- L'^iN'- + f) 

The action ([ISD rewritten in the Hamiltonian form becomes 

S= J (PlL + PrR- NH - N''Hr)drdt + Jtt^- 

I+II ^ ' 



with 



H 




(18) 



where N, , N and are Lagrange multipliers in the Hamiltonian formal- 
ism. The system of constraints contain two surface constraints in addition 
to usual Hamiltonian and momentum constraints of the ADM formalism. 
ADM constraints: 

{ ll°o (-) 

Hr = n + L[PL] = 

^ ^ R[R] , /^^2 I ^2/r2_n ^ ^ 



Shell constraints: 
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4 Kuchar variables. 



In the paper 0] Kuchar proposed canonical transformation of the variables 
{R, Pr, L, Pl) to new canonical set {R, Pr, M, Pm) in which Hamiltonian and 
momentum constraints given by (0) are equivalent to the very simple set of 
constraints : 

M' = 

The idea is to use the Schwarzschild anzatz for the space-time metric 
instead of (|l]): 



(21) 



ds^ = -F{R, m)dT^ + 



F{R,m 



(22) 



where T,R and m are some functions of (r, t) and F{R,m) = 1 — 2Gm/ R, 
and m, in general, is a function of r, m=m(r). Equating the two forms of 
the metric (|l|) and (p2D we obtain the transformation between the two sets 
of dynamical variables. 

The explicit form of the transformation is 



Pl 
R 

Pr 



(R 



l\2 



pp2 

m 



y\2 



F 

RFPm 

R 

p , Pm , FPm , {RFPrnYRR' -RFP^RR 



grfUr'Y / f-fp^ 



(23) 



where P„ 



-r. 



The Liouville form 



e 



PrR + PlL 



can be expressed in the new variables as follows: 



(24) 



e 



d 



1 



Pmfn + PrR + — \LPl + -^RR' In 



+ 



I 2G 



RRln 



dt 

RR' + LPlG 
RR! - LPlG 



2G 



RR' — LPj^G 



RR' + LPlG 



(25) 
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When there is no shell the total derivatives in ( pSj) give rise to some surface 
terms at infinities. As shown by Kuchar the appropriate falloff condi- 
tions at infinities make the last surface term from zero. Then it follows 
from ( pSD that {R, Pr, m, P^) form a canonical set of variables and equa- 
tion (|23|) describes a canonical transformation between {R, Pfi, L, P^) and 
{R,PR,m,Pm) . 
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5 Canonical variables and Hamiltonian con- 
straint on the shell. 



5.1 Shell variables. 



In the presence of the thin shell the situation is different. Surface terms now 
should not be neglected. 



Let us do the transformation (^31) in regions I and II of our space-time. The 
Liouville form of our Hamiltonian system has the form 



e = J PrR + PlL+ I vrr. 



After integration the total derivatives in 
Liouville form on S: 



e 



PrR + P^fn + 



LPr + —RR! In 



(26) 

give some contribution to the 
RR! — LPj^G 



RR' + LPlG 



RR' - LPlG 



2G 

+ J nr 
s 



RR' + LPlG 



fdt 



P^m + PrR + J pf + J Pj^R 
i+ii s s 



(27) 



where we denoted 



P 



R 



^ + L [Pl] 



20 



Rln 



RR' - GLPl 
RR' + GLPl 



(28) 



and made use of the identity 
d 



R 



dt 



R{t,r{t)) = {Rit,r) + R'{t,r)f{t)) 1=^ 



(29) 



We see that this canonical transformation involves all the set of coordinates in 
the phase space 11 = {{R{r,t) , PR{r,t) , L{r,t) , PL{r,t) , f(t) , fr{t))} according 
to the formulae (|23|) and (pSf ). Moreover it introduces additional pair of 
canonically conjugate variables {R, Pj^) on the shell. 
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In both inner and outer regions I and II constraints are simplified due to 
the canonical transformation as it was in the absence of the shell (0). The 
surface momentum constraint Hr = (^) takes the form 

p = (30) 
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5.2 Shell constraint. Special case. 

In order to investigate the form of Hamiltonian surface constraint in new 
variables let us choose coordinates (r, t) in a specific way. If coordinate lines 
t = const cross the shell perpendicularly then 



f = -A^*" 
This means according to (|T^) that 

TT = 

Then the surface momentum constraint gives 

GL [Pl] = R[FPm] = 

Let us denote 



7= FPm/L 



We have from 



R 



2G 



R 



In 



R' 



7 



R' + l 



R 
G 



In 



Kut + 1 



out 



From 



we get 



L 



— = aJ F + •y 



(31) 

(32) 

(33) 
(34) 

(35) 

(36) 



where cr = ±1 is the sign function taking its values according to whether 
radii increase in the outward normal direction to the shell ( a = +1 ) or 
they decrease ( ct = — 1 ). Using (^) and (|36|) we could find the expression 
for 7. Then from ( p6D we could determine the jump of R' across the shell 
surface as a function of P^, R, rriin and rriout- Substituting this into the surface 
Hamiltonian constraint (pOf ) we obtain the following expression 



H 



(JinR 



Pout 



Fin exp ' ^ 



Rj\ 



out 



Fin exp 



R 



-M = 
(37) 
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5.3 Shell constraint. General case. 

In the previous subsection we managed to derive the shell constraint by 
finding so some continuous variable, namely 7 (Eqn.^lD.Now we will try to 
do the same trick in general case vr 7^ 0. Let's consider the full time derivative 
of the shell radius. 

R= ^R{r{t) ,t) = R + R'r (38) 
Using the definition Pl = — ^ (-R — R'N^^ we get: 

R=-^^^pi + R'(r + N'-) (39) 
R 



Remembering that 



we can find 



ML^ (A^^ + r) 

IT 



nN 



Eqn.( R9p now reads 



It is now easy to see that the jump of R across the shell is a linear combination 
of constraints 

where 

TT 

X = — ^= 



To go further we need Eqn.(??) which we now rewrite as follows 

gPr (l + Gf (1 + x), ^ , y. n+y) 

l + Gf (1 + x)) « + i/„.*(l + x) 

/ out 
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where a = ^ and y = G^. 

The next step is to find the relation between a and R. From the definitions 
of Pl and R we have 

R'^ ^ o a u 
— = F + y\ = - + ^ 

XX 

. Solving this for y we get 

y = — ^ (43) 

Substituting this expression back into Eqn. ^ we obtain after simple 
algebraic transformations 



= ' ^^-^f+J} ^ (44) 

Z - (Touty/Z^ + Font 

where z' = ^ (45) 

1 - X 

Making use of the momentum constraint we can wright the jump of y as 
follows 



cToutXya' + (1 - X') Font = (^inXya' + (1 - X') Fin - G-^xVl - X' 

(46) 

Note that this is just the Eqn. (|^ if we choose proper time (substituting 
R for z'.) 

Squaring this equation we get 

a- + F- = + ^ 



(^outVz'^ + F, 



o«* 2MG 2R 



(47) 



z = ±. 



\\2Mg) -2^Fout + F,r.) + ^^ (48) 
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Finally we get for the shell constraint 

R [F] MG 



^ R\F\ MG ^ > 



2MG 2R 

This constraint can be rewritten in a more elegant way 



IVPG^ _ M _|_ MG 
4R2 AI 2R 



e 2R — 



where the Schwarzschild anzatz was substituted for F's . 



(50) 



It can be easily shown that the above constraint is equivalent to that one 
derived in the previous subsection. And this proves that the Hamiltonian 
constraint ( P7| ) is valid for nonzero values of vr as well. We see that the 
only remained classical constraint of the shell Hamiltonian dynamics can be 
written in various equivalent forms (e.g. , Eqns (|37|) , (^) or (^Of). But, of 
course, quantum mechanically all these forms are no more equivalent. So , 
we need some criteria to choose among them. One of this criteria well be 
considered in the last section. In what follows we will consider the following 
squared version of the Hamiltonian constraint (|37D as the suitable classical 
counterpart for the quantum constraint for the wave function 

GPr GPA M^G^ 



C = Font + Fin - \l Fout\lFin (^cxp + exp ^ j - ^'^ ^ (51) 

The Hamiltonian constraint (^) was derived under the assumption that 
both Fin and Fout are positive . It is possible, of course to derive analogous 
constraints in T-j--regions, where F < 0. But, instead, we make the following 
substitution 

and consider this function as a function of complex variable. Then the point 
of the horizon F = becomes a branching point , and we need the rules of 
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the bypass. We assume the following 



_pi/2 ^ \F\eJ-'^ 
(f) = in i?_|_-region 
= 7r/2 in T_-region (53) 
(j) = IT in i?_-region 
(/) — — 7r/2 in T+-region 



for the black hole case , and 



(f) — TT in i?+-region 

= — 7r/2 in T_ -region 

= in i?_-region 

= 7r/2 in T+-region 



(54) 



for the wormhole case. The reason for such analytical continuation is that 
we are able to get the single equation on the wave function \E' which cov- 
ers all four patches of the complete Penrose diagram for the Schwarzschild 
spacetime.Some consequences of this fact will become evident in section 6. 
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6 Quantized spherical gravity with thin shell. 



We now turn to the Dirac constraint quantization procedure . 

It is convenient to make a canonical transformation from {R, P^) to 

Ps = Rl^ 

where Rg is the gravitational radius of the shell. Dimensionless variable 5* 
is the surface area of the shell measured in the units of horizon area of the 
shell of mass M. 

The phase space of our model consist of coordinates {R{r) , Pji{r) , m{r) , 
Pm{r) , S , Ps , r , pr) {—oo,r — e)[j{r + e,oo). Then the wave function in 
coordinate representation depends on configuration space coordinates: 

^ = ^{R{r),m{r),S,r) (56) 

and all the momenta become operators of the form 

pRir) = -I 5/ 6R{r) Pm{r) = -i 5/ 5m{r) 

Ps = -I d/ ds Pr = ~id/df ^ ^ 



ADM and shell constraints (|T9D and ([20| ) become operator equations on \Ef. 
The set of ADM constraints is equivalent to the set of constraints (^Tj) in 
Kuchar variables which could be easily solved on quantum level. Indeed, in 
the regions I and II the equations 

r d^/dRir) = 

\ M'(r)^ = ^ ^ 

express the fact that wave function does not depend on R{r) and the depen- 
dence on M(r) is reduced in each region I and II to = S{M — M±) where 
M± defined in the regions I (-) and II (+) do not depend on r. They equal 
to Schwarzschild masses in the inner and outer regions Min and Mout in (|50|) . 

The set of shell constraints (pO]) impose further restrictions on First 
of them takes the form 

d^/ df = (59) 
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in new variables according to (pO|). So the only nontrivial equation is the 
shell constraint ( PD (or (^Of) , (|50|) which are classically equivalent to the 
Eqn.(0)) 

(7 (m+, m_, 5, -ihd/dS) = (60) 
vl> = vl>(m+,m_, S) 

The operator C contains the exponent of the of the momentum Ps- This 
exponent becomes an operator of finite displacement when Ps becomes dif- 
ferential operator: 

GPfl Ps _.™p'_e_ 

e H = e 2GM^ \E' = e *"M^as\l> = vl>(77i_,_, m_, S" — ^i) (61) 

where mp; is Plank mass and ^ = |(-^)^ 

The constraint C is nonlinear in 5* and Ps so the question of ordering 
should be solved when replacing dynamical variables by operators. It is 
proposed in 0to choose the symmetric ordering 

A{S)B{Ps)^ l/2{A{S)B{-tnd /dS) + B{ind /dS)'A{S)} (62) 

where A and B are some functions of S and Ps respectively and A denotes 
complex conjugation. 

With operator ordering (|^) we must add to the constraint (^) the com- 
plex conjugate part . 

The constraint C becomes an equation in finite differences if we express 
R through 5* and substitute the expression (??) to the differential operator. 
Finally we get 



We have mentioned already that the classically equivalent constraints give 
inequivalent quantum theories. This is well known fact. We suggest that the 
criterion to choose the correct quantum theory is the behavior of the wave 
functions in the quasiclassical regime. In our case this means the large black 
holes limit. Indeed , the parameter = |("~f")^ becomes small for large 
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and the expansion with respect to this parameter is equivalent to 



the expansion in Planckian constant h {nipi = ^h/Gc). In the next Section 
we will consider this quasiclassical limit and show that our choice for the 
quantum constraint is a good one (at least in the case of one thin shell). At 
the and of this Section we would like to make an important remark. Our 
quantum equation ^ (which is just a Schroedinger equation ) is the equa- 
tion in finite differences rather than differential equation, and the shift in 
argument is along an imaginary axis. In the case of differential equation we 
we require the solution to be differentiable sufficiently many times . Simi- 



larly, we have to demand the solutions of our finite differences equation |63 



to be analytical functions. This condition is very restrictive but unavoid- 
able. Our previous experience (see 0) shows that it is the analyticity of the 
wave functions and not the boundary conditions that lead to the existence of 
the discrete mass (energy) spectrum for bound states. How it works in the 
quasiclassical regime we will see in the next Section. 
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7 Large black holes. 



The finite difference equation ( |63D becomes an ordinary differential equation 
in quasiclassical limit which is the same as the limit of large ( m ^ rupi ) 
black holes. Indeed the parameter of finite displacement of the argument of 
\l/ in (^) C, = {fnpi I M) becomes small and we could cut the Tailor expansion 

^{S + ^z) = ^{S) + z^^'iS) - ^^"{S) + ... (64) 

at the second order. 

Now we will analyze the behavior of the solutions of equation ( |63D in this 
quasiclassical limit at singular points. 

We will restrict the consideration to the case of flat inner region m_ = 0, 
so we denote m+ = m. It is convenient to redefine the area variable S so 
that 

' = (2^ = :rI (65) 



^9 2R 

area is now measures in the units of horizon area of a black hole with 
Schwarzschild mass M {Rg is its gravitational radius). In these units the 
displacement parameter 

C = ^ (66) 

m 

and the equation (^) reads as 

e^^^l {^{s + <) + $(s - <)) + e-**y^(s + <)^(s + <)+ .g^. 
e-<^y^(s - <)vl/(s - <) = 2(F - ^)v|/(s) 

where is the phase of F^/^ . It should be chosen in different R- and T- 
regions according to the arguments of section 5. In the last formula we must 
take the Tailor expansion on ^ up to second order. 

v^l^^^, ^ ^{S)±^'{S)Cz-^^"{S)... 



S±Ci 



(68) 
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This leads to ordinary differential equations of second order, which are differ- 
ent in _R+, R_,T+ and T_ regions due to the different values of the phases in 
Eqn.(^). The interesting for us singular points of these differential equations 
are 

S = oo and S = 1. (69) 

In the quasiclassical limit our requirement of the analyticity of the solutions 
to the exact equation (^) transforms into the requirement that the branching 
points of the leading terms in the solutions to the approximate equations 
should be of the same kind . Thus, we need to keep only those terms in the 
corresponding equations that give us these leading terms. Below we consider 
the black hole case only. The results are easily translated to the wormhole 
case. 

The singular point 5 = cxd in the region i?+ lies in a classically forbidden 
region as far as we restrict ourselves with bound motions of the shell only. 
In order to analyze the behavior of \E' in this region we should take (|6^) with 
= and expand all the quantities in terms of y, where s = (1 + y)"^ The 
result is 

1 1 / 1 MM 

%y--% + T2{l- — + -{2-—)]^ = 70 
y ( \ m 2y m J 

The leading term of the solution is 

1 — 2 ~ 2 
_ _ rn 

^ ^ y2 A^iC exp(-/iy), (71) 



1 /m^ 



For another singular point in i?+ region , that is for 5 — » 1 + we have 
v^..-3.v^. + ^(l-^U = (72) 



with leading term 



y = v^, s > C, y > C, C < 1 
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Comparing the types of the branching points at s — > oo and s — > 1 + 
we can conclude that 



2- — 

m 



— n, n — integer (74) 



V m 

This is the first quantization condition. We will not consider here the worm- 
hole case . Note only that , as can be shown , the positive values of quantum 
number n correspond to black holes while negative n correspond to worm- 
holes. 

In the T_-rcgion (which classically is a region of inevitable contraction) , 
i.e., for s ^ 1 — (s = (1 + y)"^, y < 0, C ^ 1^1 ^ 1) have first order differ- 
ential equation (due to complex conjugation introduced earlier the leading 
terms containing second derivatives of the wave function cancel each other) 
with following leading term in the solution 

which is just the ingoing wave as it should be expected for the quasiclassical 
limit in the region of the inevitable contraction. That is why we have chosen 
the of the function (F)-*^/^ = e*'^|Fp/^ in the TLregion (Which classically the 
region of inevitable expansion) the choice oi (f) — — | leads to the outgoing 
wave as a solution. Note also that to the our requirement of the analyticity 
the solution (leading term) in the T_rcgion should be the analytical contin- 
uation of the solution in the i^+region. And we see that this is indeed the 
case. 

We do not consider here separately the asymptotics in i?_ region near 
the horizon (s — > 1 -|- 0) because it differs from the corresponding solution in 
i?+-region only by the sign in front of the second term. 

Let us now turn to the asymptotics of the solutions in i?_-region for 
s — > oo. Due to the minus sign in front of F^/^ the equation for the wave 
function in a i?_-region is quite different from that in a i?+-region 

1 1 / \ 
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The leading term of the asymptotic is now the following 



y 



exp-(-y 



(77) 



Note that the falloff in the i?_-region is much faster than it is in the 
region. This is a quite reasonable result because it means that the quantum 
shell in the black hole case can penetrate into the i?_-region (which is com- 
pletely forbidden for the classical motion) but the probability of such an 
event is negligible small. 

And , again, comparing the types of the branching points at s ^ 1 + 
and s — »• oo in the i?_-region we get 



m 



- 1 



+ p, p = positive integer 



(78) 



8C 2 

The appearance of the4 second quantum number is rather surprising but it 
allows some explanation. We discuss this point in the last Section. 
Combining (73) and (|78|) we get 



. 2 

m 



1 



>3/2 



l + 2p 



m 



n 



(79) 



and 



m 



-m 



m 



pi 



(80) 



For p 3> \n\ we have 

m ^ 2y/p rUpi (81) 

This corresponds to the shells of large mass whose mean value radius is rather 
close to the horizon. In the opposite case , p <^ n, 



m 



V2{1 + 2pf'^n^l^mpi 



(82) 



26 



This corresponds to the massive shells with mean radius very far from hori- 
zon. 

This the end of this Section we would like to consider the behavior of 
the solutions in the vicinity of the horizon (sub-Planckian deviation ), where 
\y\ ^ C (-s ~ !)• To be specific we will be interested in the solutions -R+ and 
T_ regions. The expansion ( |68D is no more valid for the function F^/'^{s±iC,) 
but it is still valid for the wave Keeping the leading terms only we have 
now 



with the solution 




1 - 



4m^ 



^(s) = 



\& ~ e 



C \ qC 



-5/2 



(1 



Aw? 



(83) 



(84) 



The coefficient a equals to 1 in the i?+-region and to imaginary unit i in the 
T_-region. 

In the i?+ region 



±i 



4 



Am' 



(85) 



and we have superposition of two waves (ingoing and out outgoing) with 
relatively equal amplitudes, 
in the T_-region 



(1--^)(1 + ^) 




(86) 



The existence of two waves in the T_-region reflects the quantum trembling 
of the horizon . But the outgoing wave is enormously damped relative to the 
ingoing wave (of course , in the T_- region the situation is exactly inverse 
one). It is this damping that cause (in a quasiclassical regime) existence of the 
single ingoing wave in the T_-region at the distances larger than Planckian. 
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8 Discussion 



In the concluding remarks we would like to discuss the obtained results. 

(1) We constructed the geometrodynamics for the spherical gravity with 
self-gravitating thin dust shell as a source. Such a shell provides us with the 
only dynamical degree of freedom which is otherwise absent in pure spherical 
gravity. We managed to separate the field canonical variables describing the 
gravitational field outside the shell from the single pair of the shell canon- 
ical variables which are just the radius of the shell and the corresponding 
conjugate momentum. The Hamiltonian constraint is derived for the shell 
which depends only on the invariants of the inner and outer parts of the 
manifold and on the parameters of the shell. The quantum functional is 
subject to both quantum ADM (field) constraints and the shell constraint. 
After solving the (trivial) quantum field constraints we are left with the 
functional which depends on the inner and outer masses of the correspond- 
ing Schwarzschild manifolds and is function of the radius of the shell. Thus , 
the functional becomes a wave function , and the remaining shell constraint 
is just a Schroedinger equation for this wave function. 

(2) The obtained Schroedinger equation is the equation is the equation 
in the finite differences rather than differential equation. And the shift is 
along an imaginary axis . Dealing with the differential equations we always 
require (or assume) that the solutions should be sufficiently differentiable. 
Analogously in our situation we must require solutions to be analytic func- 
tions except some isolated points . Our equation has branching points , so , 
the solutions will have branching points as well. 

(3) The Schroedinger equation we obtained contains , as a coefficient 
function , the square root of the F = 1 — which is invariant function 
of the Schwarzschild solution. This function is positive outside the event 
horizons on the both sides of the Einstein-Rosen bridge and -R_- regions) 
and it is negative beyond the horizons in the T+-region of the inevitable 
expansion and in the T_-region of the inevitable contraction. We suggest to 
consider the square root as a function of complex variable ,(-F)^^^ acquiring 
different phases in patches of the complete Schwarzschild manifold. The 
aim of such a procedure is twofold. First, it allows us to obtain a common 
wave function covering the hole Penrose diagram (R-^ and T-t-regions). And 
, second, we remove the double cover degeneracy when the same value of 
radius R corresponds to two different points , on in the i?+(T+) -region and 
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the other is in the i?„(T_) -region. Now this different points with the same 
value of radius he in different sheets of the Riemannian surface. 

(4) The requirement of the analyticity of the wave function on the cor- 
responding Riemannian surface means that the branching points should be 
of the same kind in order to be connected by cuts. In another words , the 
number of the Riemannian sheets emerging at the branching points must be 
the same. 

In our case we have the branching points at infinity and at the horizon. 

But , we have two different horizons. One of them separating, say and 
T_-regions , lies on the one sheet and the other one , separating T_ and 
regions lies on another sheet of the Riemannian surface . Thus comparing 
the branching points at infinity and at the horizon in i?_|_-region we obtain 
the first quantum number characterizing the mass spectrum of the system 
in question and, comparing these points in the TL region we obtain the 
second quantum number. Thus , the mass spectrum of the black holes and 
wormholes should depend on two quantum numbers . (Note that using the 
same method one can obtain the famous spectra like oscillator, hydrogen 
atom and so on). 

(5) It is well known that the classical theory may give rise to different 

inequivalent quantum theories. The origin of this phenomenon is a non- 
commutativity of dynamical variables and their conjugate momenta. The 
investigation of the quasiclassical limit helps to choose the "correct" quantum 
version . In our case the quasiclassical regime coincide with the limit of 
large (comparing to the Planckian mass) black holes. The finite differences 
equation can now be expanded in series with respect to the small parameter 
and we can cut the series to obtain the differential equation. We showed that 
our choice of the quantum Hamiltonian give a good quasi- classics (ingoing 
wave in the TLregion , and outgoing wave in the T^-region, as should be 
expected). Moreover we showed that the black hole and wormhole discrete 
mass spectra is determined by two quantum number making this spectra 
quasicontinuous. This resembles the appearing of the fine structure due to 
removing some degeneracy (in our case it is a double covering degeneracy). 
In the ordinary quantum mechanics we are used to the fact that the number 
of quantum numbers equals to the numbers of degrees of freedom. Prom 
the first sight we have only one dynamical degree of freedom in our problem. 
Indeed, the motion of the spherically symmetric thin shell is described by the 
radius as a function of time in Lagrangian picture and by the radius and its 
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conjugated momenta in the Hamiltonian picture. But actually we have two 
different parts of the Penrose diagram with the same value of radius , namely 
and i?_-regions (and T+ and T_-regions). And by our consideration they 
lie on different sheets of the Riemannian surface. Of course, classical motion 
is forbidden in i?_region in the black hole case (in i?_|_region in the wormhole 
case). But in the quantum theory such motion is allowed. Thus , we have 
in fact two degrees of freedom . And this is just the origin if the second 
quantum number. 

(6)This last remark concerns the problem of small black holes with mass 
about Planckian mass or smaller. The shift in our equation in finite differ- 
ences is of order of the horizon size or even larger. This means that for small 
masses our equation does not feel the very existence of the horizons. And 
it gives some hope that there are no black holes at all with masses smaller 
than the Planckian mass. 
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q„ =+1 



Fig. 2. Black hole case. 



M>0 




Fig. 1. Penrose diagram for Schwarzchild black hole 
Dashed lines are curves of constant radius. 



»i,„ > m„„, M>0 



Fig. 3. Black hole case. In- and out- regions are interchanged. 



q„ = -1 



Fig. 5. Unphysical. 



M<0 




Fig. 4. Wormhole case 



